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Abstract. The reorientation of one small paramagnetic molecule ( spin probe ) in 
glassy polystyrene ( PS ) is studied by high-field Electron Spin Resonance spectroscopy 
at two different Larmor frequencies ( 190 and 285 GHz ). The exponential distribution 
of the energy-barriers for the rotational motion of the spin probe is unambigously 
evidenced at both 240 K and 270 K. The same shape for the distribution of the energy- 
barriers of PS was evidenced by the master curves provided by previous mechanical 
and light scattering studies. The breadth of the energy-barriers distribution of the 
spin probe is in the range of the estimates of the breadth of the PS energy-barriers 
distribution. The evidence that the deep structure of the energy landscape 
of PS exhibits the exponential shape of the energy-barriers distribution 
agrees with results from extreme-value statistics and the trap model by 
Bouchaud and coworkers . 
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1. Introduction 

The study of glassy solid dynamics is a very active one PIE] • Here one is interested in the 
temperature range which is, on the one hand, well below the glass transition temperature 
T g to neglect aging effect and consider the glassy system as one with constant structure 
and, on the other hand, high enough to neglect tunneling effects governing the low- 
temperature anomalies of glasses. In this regime the dynamics is thermally activated in 
the substructures of the minima of the energy landscape accounting for various subtle 
degrees of freedom [3] . It has been shown that in the glass the temperature dependence 
of the energy barrier distribution g(E) is only weakly temperature-dependent [HE]- 

The shape of the energy-barriers distribution g(E) in glasses has been extensively 
investigated via experiments 011113013110111111121131111, theories [13 HQ El CHI 
[13 E3 EH E2] and simulations j2S] • Basically, two different distributions are usually 
recovered, the gaussian distribution IH 13 HU 113 [HI E3 and the exponential 
distribution [TTJ [T21 HS1 HE1 UHl UHl E01 EH E2! - The convolution of these two distributions 
[21] as well as the truncated Levy fligth, i.e. a power law with exponential cutoff, 
resembling the stretched exponential ^1] were also considered. 

It is interesting to relate g(E) with the density of states, i.e. the distribution of the 
minima of the energy landscape. On the upper part of the landscape, being explored 
at high temperatures, the Central Limit theorem suggests that the density is gaussian 
[13 E3j- At lower temperatures the state point is trapped in the deepest low-energy 
states which are expected to be exponentially distributed following general arguments 
on extreme- value statistics leading to the so-called Gumbel distribution [20]. Random 
Energy models |15 3 and numerical simulations [22] support the conclusion. The barrier 
height E during the jump from one state with energy E\ to another state with energy 
E 2 has been modelled by the linear combination E = aE 2 + (1 — a)E\ [21] • In the case 
of trap models ( a = ) the minima and the energy barriers have the same distribution. 
This also holds true for annealed disorder and a > ^3E]- Moreover, one notes that 
the linear combination of independent gaussian variables is gaussian too and the linear 
combination of independent exponential variables has exponential tail. All in all strict 
relations are expected between g(E) and the density of states. 

If the average trapping time r before to overcome the barrier of height E at 
temperature T is governed by the Arrhenius law, 



k being the Boltzmann's constant, the distribution of barrier heights induces a 
distribution of trapping times p{r) . The explicit form of p(r) for a gaussian distribution 
of barrier heights with width <je is the log-gauss distribution ( LGD ) 



r = T exp(E/kT) 



(1) 





Energy landscapes and High-field ESR 



3 



° — &E/kT is the width parameter. p{r) is maximum at tlgd- If the distribution of 
barrier heights is exponential with average height E 

g{E) = =exp(~) (3) 

p(r) is expressed by the power-law distribution ( PD ) 

/ s J if r < t pd 

Ppd{t) = < x _ {x+1) 4 

I XTp D T (X+L> it r > T PD 

with x = kT/E and Tpp, = r . If the width of the energy-barriers distribution is 
vanishingly small, a single trapping ( correlation ) time ( SCT ) is found and both Plgd 
and ppd reduce to: 

Psct(t) = 5(t - t sct ) (5) 

The shape of g(E) has been usually studied by measuring different susceptibilities 
x"( v ) arising from either collective or single-particle response. If the susceptibility 
follows by a static distribution of activated relaxation times in the presence of a wide 
distribution g(E) it follows 



X"{y) oc Tg(E), E = kTln(l/27rvT ) (6) 

The above equation shows that the shape of x"{ v ) yields the shape of g(E). The 
conversion factor between the frequency and the energy scales is the unknown attempt 
frequency l/r which is usually treated as one adjustable parameter to set both the 
width and the location of g(E). The dielectric spectroscopy provides a convenient 
frequency range to recover the full shape of x"{ v ) and then of g(E) [H IHJ- In other 
cases, e.g. light scattering [12^ and mechanical relaxation ^3], the accessible frequency 
range is limited and g(E) is recovered by building suitable master curves assuming the 
time-temperature superposition principle. Nuclear Magnetic Resonance ( NMR ) offers 
an alternative procedure to get information on g(E) by analyzing the NMR lineshape 
in terms of two weighted components [2l|2Hj. However, an adjustable conversion factor 
between the frequency and the energy scales is also needed. 

The use of suitable probes to investigate the secondary relaxations in glasses by 
NMR El El EE!, Electron Spin Resonance ( ESR ) |2H EHl W\ and Phosphorescence 
|3H] studies is well documented. In spite of that efforts, the relation between the probe 
motion and the host dynamics is usually not obvious with few exceptions which,notably, 
involved relatively small and nearly spherical probe molecules [2E]- It was also noted 
that small molecules, e.g. xanthone and benzophenone, are more sensitive to shorter 
segmental motions occurring at lower temperatures [30J. 

During the last few years continuous-wave ( CW ) and pulsed High-Field ESR 
( HF-ESR ) techniques were developed involving large polarizing magnetic fields, e.g. 
Bq = 3T corresponding to Larmor frequencies about 95GHz ( W band ), [3T1 IH21 EH| • 
HF-ESR is widely used in solid-state physics jHlj, biology [SHI ESI EH ESI EH] and polymer 
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science |4U[ I4*T] W2\ I43j . One major feature is the remarkable orientation resolution [4*3] 
due to increased magnitude of the anisotropic Zeeman interaction leading to a wider 
distribution of resonance frequencies [313 • 

In the present paper one demonstrates that the rotational motion of suitable small 
guest molecules ( spin probes ), as detected by HF-ESR , is an effective probe sensing 
the energy-barrier distribution g(E) of glassy polystyrene ( PS ). The choice of PS 
was motivated by a number of studies of g(E) by Raman ^1] and light scattering 
and mechanical relaxation [THj which evidenced its exponential form ( eqEl ). The 
distribution of barriers to be surmounted by suitable probes in PS was also studied by 
NMR [23]. The secondary relaxations of PS were investigated by phosphorescent probes 



The paper is organized as follows. In Sec. [2] the background on ESR is presented. 
Experimental details are given in SecJH] and the results are discussed in Sec. [31 The 
conclusions are summarized in Sec. EJ 

2. ESR Background 

The main broadening mechanism of the ESR lineshape is the coupling between the 
reorientation of the spin probe and the relaxation of the electron magnetisation M 
via the anisotropy of the Zeeman and the hyperfme magnetic interactions. When the 
molecule rotates, the coupling gives rise to fluctuating magnetic fields acting on the spin 
system. The resulting phase shifts and transitions relax the magnetisation and broaden 
the resonance [3HI15]. Because of the roughness of the energy landscape and the highly 
branched character of the free volume distribution, one expects that the small spin probe 
undergoes jump dynamics in glasses . One efficient numerical approach to calculate the 
ESR line shape in the presence of rotational jumps is detailed elsewhere |41j . 

The rotational correlation time r c , i.e. the area below the correlation function of 
the spherical harmonic Y2,o, is [27, 



where and r* are the size of the angular jump and the mean residence time before 
a jump takes place, respectively. In the limit (f> « 1 the jump model reduces to the 
isotropic diffusion model with r c = 1/6D = t* / cf) 2 where D is the rotational diffusion 
coefficient. 

The occurrence of a static distribution of correlation times in glasses suggests to 
evaluate the ESR line shape L(B ), which is usually detected by sweeping the static 
magnetic field B and diplaying the first derivative, as a weighted superposition of 
different contributions: 



EH- 



T* 



(7) 



1 - 





(8) 
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L(B , t c ) is the EPR line shape of the spin probes with correlation time r c and p(r c ) is 
the t c distribution. 

It must be pointed out that in the presence of wide distribution of correlation times 
the ESR lineshape ( eq. ) cannot be simplified as a sum of few leading components. 
This is rather different with respect to NMR where 'two-phase spectra' stemming from 
'fast' and 'slow' molecules are observed [2J. If, on one hand, that feature leads to more 
elaborated numerical work, on the other hand it ensures that the ESR lineshape exhibits 
a much richer variety of details being markedly affected by the spin-probe dynamics 
which provides helpful constraints to the fit procedures and results in an improved 
information content of the spectra. 

3. Experimental 

PS was obtained from Aldrich and used as received. The weight-average molecular 
weight is M w = 230kg jmol and T g = 367 K. The free radical used as spin probe 
was 2,2,6,6-Tetramethyl-l-piperidinyloxy (TEMPO) from Aldrich. TEMPO has one 
unpaired electron spin S — 1/2 subject to hyperfine interaction with the nitrogen 
nucleus with spin 1=1. The sample was prepared by the solution method [IB] by 
dissolving TEMPO and PS in chloroform. The solution was transferred on the surface 
of one glass slide and heated at T g + 10K for 24h. The spin probe concentration was 
less than 0.08% in weight. Appreciable broadening of the ESR line shape due to the 
spin-spin interaction are observed for concentration larger than 0.2% in weight. No 
segregation of the spin probe was evidenced. Samples aged at room temperature for six 
months exhibited no appreciable changes of the ESR lineshape. 

The ESR experiments were carried out on the ultrawide-band ESR spectrometer 
which is detailed elsewhere In the present study the spectrometer was operated at 
two different frequencies, i.e. at 190GHz and 285GHz. The multi-frequency approach 
ensures better accuracy to determine the spin-probe dynamics. The sample of about 
0.8cm 3 was placed in a Teflon sample holder in a single-pass probe cell. All spectra were 
recorded and stored in a computer for off-line analysis. 

4. Results and Discussion 

The numerical simulation of the ESR lineshape much relies on the accurate 
determination of the magnetic parameters of the spin probe. To this aim, one profited 
from the enhanced resolution of HF-ESR and measured the ESR lineshape of TEMPO 
in PS at low temperature where the rotational motion of the spin probe is very slow 
and the resulting lineshape approaches the so called 'powder' or 'rigid' limit lineshape 
|45j . The results at 50K are shown in Fig^ It is seen that the lineshapes at both the 
operating frequencies are well fitted by the SCT model with a single set of magnetic 
parameters. The small discrepancy between the simulation and the lineshape at low 
magnetic field was already noted in other studies [32] ■ Remarkably, no distribution 
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of correlation times is evidenced at 50if. At this stage a complete discussion of this 
finding needs the detailed investigation of the temperature dependence of the rotational 
dynamics of TEMPO in PS. This is beyond the purpose of the present Letter and will 
be discussed thoroughly in a forthcoming paper [38] . 

FigE] compares the ESR lineshapes at IWGHz of TEMPO in PS at T = 270K = 
T g — 97K with the SCT ( left, two adjustable parameters t SC t,4> ) an d LGD models 
( right, three adjustable parameters r LGD ,(j) and a ). For a given jump angle <p the 
remaining one ( SCT ) or two ( LGD ) parameters have been adjusted. It is apparent 
that the SCT model is inadequate. Large disagreements were also found when fitting 
the ESR lineshape at 285GHz ( not shown). FigE] shows that the LGD model yields 
better agreement for large jump angles. However, closer inspection reveals deviations 
in the wings of the ESR lineshapes at both low- and high-field. The ESR absorption 
of these regions is mainly contributed by spin probes with their x molecular axis (low- 
field) or z axis ( high field) being parallel to the static magnetic field 03] • Being the 
magnetic parameters of the spin probe precisely measured, it may be shown that the 
above disagreements are clear signatures of the overestimate of the jump angle size 
[48] . The best-fit value of the width parameter a = 1 of the LGD distribution (eq. EJ) 
corresponds to the width a^/k = 270K of the energy-barrier distribution of TEMPO 
in PS. Rossler and coworkers found by NMR cr^/fc = 276K for hexamethylbenzene 
in PS in the temperature range roughly 150 — 300K [21]. Even if TEMPO and 
hexamethylbenzene have similar sizes their shape is rather different, then it is quite 
reassuring to note that the rotational motion of different small probes investigated by 
different techniques lead to comparable results. 

Fig 0] compares the PD model ( eq. HJ, three adjustable parameters tpd,4> and x ) 
with the ESR lineshape at both 190 and 285GHz of TEMPO in PS at 270K. Again, 
having fixed the jump angle, one adjusted the other two parameters. It is clearly seen 
that the agreement is quite improved, especially for small jump angles. Fig. HJ provides 
the same comparison at 240K for the case <fi = 20° which results in the best agreement. 
From FigsJUandHJit is seen that at both 2A0K and 270K the best-fit value of 
r PD is different for the two operating frequencies 190GHz and 285GHz with 
tpd(190GHz) < Tprj(285GHz). The uncertainty on tpd is roughly less than 20%. 
Constraining Tpp> to having the same value at 190GHz and 285GHzin the fit 
procedure results in larger deviations between the experimental lineshape 
and the numerical results. The decrease of Tpp> with the operating ESR 
frequency has been noted by us also at frequencies lower than 190GHz and 
will be discussed elsewhere |48j. Without going into details, it has to be 
ascribed to the relation between the Larmor period ( roughly the inverse of 
the operating frequency ) and the sensitivity of the ESR spectroscopy to the 
fast rotational dynamics. 

Table [0 summarizes the best-fit results of the PD model ( <fi = 20° ) for both 190 
and 285GHz. The width of the barrier-height distribution of TEMPO in PS is in the 
range 470K < E/k < 705K. One may wonder if this range corresponds to correlation 
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times which are effectively accessed by HF-ESR in that the longest correlation time 
of TEMPO which may be measured is T max = 100ns. To this aim, one estimates the 
maximum energy barrier which TEMPO may overcome leading to appreciable motional 
narrowing effects in the lineshape as E max = kT ln(T max /T PD ) . In the temperature range 
2AQK - 270K we get E max /k = 1500 - 1800K which is fairly larger than E, see Table 
^ This confirms that the overall ESR lineshape, expressed by the superposition given 
by eqlHl is mostly due to components which are motionally narrowed to some or large 
extent by the reorientation of TEMPO. Table Q shows that the width of the barrier- 
height distribution of TEMPO in PS increases on cooling. This points to better coupling 
of TEMPO and PS at T = 240iC However, we remind that at lower temperatures the 
width decreases ( see FigJI]). This non-monotonic behavior, as well the increase of E 
with the frequency which is also apparent in Tabled will be discussed in detail elsewhere 

EH]. 

The exponential distribution ( eq. |3] ) of the barrier-heights of PS was evidenced 
by both internal friction ^3], Raman ^1] and light scattering ^2] measurements. The 
studies converted the mechanical [13 a and the optical [TTJ [T2| susceptibilities y!'{y) by 
eqlHlto get g(E) and estimated its width as E IF /k = 760±40-K" , E Raman /k = 530±60fT 
and E LS /k = 530 ± 40K, respectively. It was concluded that, in spite of the different 
estimates of E, the corresponding three patterns of g(E) are in good agreement ( see 
ref.^2] and especially fig. 4 ). Table ^ shows that the width of the barrier- height 
distribution of TEMPO is in the range of the estimates on the width of the PS 
distribution provided by the above studies. 

One final remark concerns Tp^, i.e. the time scale tq of the activated 
trapping of TEMPO, eq. [TJ The best-fit values of r are about 120 - 250ps 
( see captions of Figs. El SI)- Previous studies about the rotational dynamics 
of TEMPO in several glassy polymers were carried out by the customary X- 
band ESR spectroscopy at 9GHz and analysed in terms of the SCT model 
[HI HOI- They yielded r = 10 2 - 10 3 ps . In particular, for TEMPO in glassy 
PS it was found r = 10 2 ps [49 J. 

5. Conclusions 

The rotational motion of the guest molecule TEMPO in PS was investigated by using 
HF-ESR at two different Larmor frequencies ( 190 and 285 GHz ). The use of a single 
correlation time ( SCT model ) was found to be inadequate to fit the ESR lineshape. 
Limited improvement was reached by considering a gaussian distribution of energy- 
barriers for TEMPO ( LGD model ) even if the best-fit width compared rather well 
with the one of molecular probes with similar size in PS as measured by NMR studies 
|24j . Assuming the exponential shape of the energy-barriers distribution for the spin 
probe ( PD model ) led to much better agreement at both 240 K and 270 K. The 
same shape was also evidenced by other studies [TTJ [T2J HB] on the distribution of the 
energy-barriers of PS which reported considerably different estimates of the width E 
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( > 40% ). The width of the barrier- height distribution of TEMPO was found in the 
range of the previous estimates of E for PS. 

The evidence that the deep structure of the energy landscape of PS 
exhibits the exponential shape of the energy-barriers distribution agrees with 
results from extreme-value statistics [20J and the trap model by Bouchaud 
and coworkers [lJJJ 121] . 
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Table 1. The width of the exponential distribution of barrier- heights E of TEMPO in 
PS at T = 2A0K and 270K as provided by the best-fit of the HF-ESR lineshapes at 190 
and 285GHz by using the PD model with jump angle (j> = 20°. Previous measurements 
of the width of the distribution of barrier- heights of PS by internal friction ^3], Raman 
[TT] and light scattering [T2| yield E IF /k = 760 ± AOK , E Raman /k = 530 ± 60K and 
E LS /k = 530 ± 40K, respectively. 

T{K) f(GHz) E/k(K) 

240 190 600 ± 36 

240 285 705 ± 42 

270 190 470 ± 28 

270 285 540 ± 32 
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Figure 1. The ESR line shapes at 190GHz (left ) and 285GHz ( right ) of TEMPO in 
PS at 50K . Magnetic parameters are: g x = 2.00994 ± 3 • 1(T 5 , g y = 2.00628±3- 10" 5 , 
g z = 2.00212 ± 3 • 10~ 5 , A x (mT) = 0.62 ± 0.02, A y {mT) = 0.70 ± 0.02, A z {mT) = 
3.40 ±0.02 The superimposed dashed lines are best- fits according to the SCT model eq. 
El with tsct = 25ns ( 190GHz ) and t SC t = 19ns ( 285GHz ). Jump angle = 60°. 
The theoretical lineshapes were convoluted by a gaussian with width w = O.IStoT to 
account for the inhomogeneous broadening. 



Figure 2. Comparison of the EPR line shape at 190GHz of TEMPO in PS at 270K 
with SCT and LCD models. Magnetic parameters and gaussian width as in FigU 
Left: Best fits according to the SCT model , eq. El with different jump angles. The 
best-fit parameters are: = 5°, tsct — 4.2ns; </> = 20° , tsct — 4.16ns; <j> = 60° 
, tsct — 3.75ns; <j> — 90° , tsct — 3.16ns. Right: Best fits according to the LGD 
model, eqEl with different jump angles. The best-fit value of the width was found to 
be independent of the jump angles and set to a = 1.0. The other best-fit parameters 
are: = 5°, tlgd = 3.6ns; <fi = 20° , tlgd — 3.6ns; (p — 60° , tlgd — 1.8ns; (j) = 90° 
, t lgd = 1.6ns. 



Figure 3. Best fits of the EPR line shape at 190GiJz ( left ) and 285GHz ( right 
) of TEMPO in PS at 270i^ according to the PD model, eqU and different jump 
angles <p. Magnetic parameters and gaussian width as in Fig^ Best-fit parameters 
as follows. Left: = h°,x = 0.6, t pd = 0.24ns; (f> = 20°, x = 0.575, t pd = 0.225ns; 
4> = 60°, x = 0.6, t pd = 0.25ns; = 90°, x = 0.63, t pd = 0.25ns. Right: (j> = h°,x = 
0.52, t pd = 0.13ns; (f> = 20°, x = 0.5, t pd = 0.13ns; <j> = 60°, x = 0.5, t pd = 0.15ns; 
(f, = 90°, x = 0.55, t pd = 0.15ns. 



Figure 4. The ESR line shapes at T = 24QK and frequencies 190GHz (panel a ) and 
285GHz (panel b). The dotted superimposed curves are numerical simulations by using 
the PD model ( eq|U) with x = 0A,t pp > = 0.25ns (panel a ); x = 0.34, t pp > = 0.12ns 
(panel b ). In both cases the best-fit value of the jump angle is <fi — 20° . Magnetic 
parameters and gaussian width as in FigE] 
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